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ABSTRACT 

We compare the results of the mark correlation analysis of galaxies in a sample from the Sloan Digital Sky 
Survey and from two galaxy catalogs obtained by semi-analytical galaxy formation models implemented on the 
Millennium Simulation. We use the MOPED method to retrieve the star formation history of observed galaxies 
and use star formation parameters as weights to the mark correlations. We find an excellent match between 
models and observations when the mark correlations use stellar mass and luminosity as weights. The most 
remarkable result is related to the mark correlations associated to the evolution of mass assembly through star 
formation in galaxies, where we find that semi-analytical models are able to reproduce the main trends seen 
in the observational data. In addition, we find a good agreement between the redshift evolution of the mean 
total mass formed by star formation predicted by the models and that measured by MOPED. Our results show 
that close galaxy pairs today formed more stellar mass ^ 10 Gyr ago than the average, while more recently 
this trend is the opposite, with close pairs showing low levels of star formation activity. We also show a strong 
correlation in simulations between the shape and time evolution of the star formation marks and the number 
of major mergers experienced by galaxies, which drive the environmental dependence in galaxy formation by 
regulating the star formation process. 

Subject headings: galaxies: stellar content — large-scale structure of universe 



1. INTRODUCTION 

Tools to extract physical information from galaxy spectra 
have evolved spectacularly in the past five years. It is now 
possible to extract accurately the past star formation (SF) 
and metallicity histories of galaxies from thei r integrated stel- 
lar light (e.g . Heavens et al. 2004; Cid Ferna ndes et al.ll2b05l : 
iPanter et al.li2007c iCid Fe rnandes e t al.ii2007l) . Further, with 
the advent of mega-spectroscopic surveys, it is now possible 
to have samples large enough to study the physical param- 
eters that drive the formation and evolution of galaxies as a 
function of galaxy properties. One obvious parameter is the 
galaxy mass, and it seems now firmly established that mass 
has an influence on the star formation rates (SFRs) of galax- 
ies (e.g. lCowie et alJil996 ). However, and open question is if 
other parameters besides mass influe nce the shape and fo rma- 
tion of galaxies (e.^ao et al. 2005 UAngulo et al.ll2007h . 

In a recent study, S heth et al.i (i200 ot used mark correlations 
of the MOPED SF and metallicities of SDSS galaxies to de- 
termine how these quantities were correlated. They found that 
close pairs ( < 1 Mpc) at z ~ O.l dominated the SF ac- 
tivity in the past relative to the mean, while today the SF is 
dominated by pairs that are separated by more than 10 Mpc. 
Similar trends were found for stellar metallicity. While this 
result fits well within the framework of galaxy "downsizing" 
it is unclear how current models of galaxy formation are able, 
or not, to reproduce it. In fact, a halo occupati on distribu- 
tion (HOD) study shows that the trends found in I Sheth et al] 
(|2006) cannot be reproduced by changing the free parameters 
in the HOD model, thus pointing toward a need for some en- 
vironmental dependence, which the current HOD paradigm 
does not incorporate. Thus, it is useful to look into numerical 
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simulations how environment can influence galaxy evolution. 

In this Letter, we explore how semi-analytical models 
(S AMs) of galaxy formati on perform at repr oducing the mark 
correlations computed in Sheth et aTl (1200 6'). We use the Mil- 
lennium Simulation and two semi-analytical recipes imple- 
mented on it to follow galaxy evolution. We find that models 
are successful at reproducing the trends observed in the SDSS 
mark study. The models therefore do well at reproducing the 
environmental trends of SF suggesting that mass is not the 
only driver of galaxy formation and evolution. 

2. STAR FORMATION HISTORY OF GALAXIES 

The observational galaxy sample used in this work was 
drawn from the SDSS Data Release 3 (DR3; Abazaiian et aD 
l2005h . We built a volume limited sample with 7'-band absolute 
magnitude Mr < -20 and redshift range 0.02 < z < 0.071, 
with a conservative limit in the apparent magnitude of < 
17.5. This samp le corresponds to the faint catalog used by 
IShefli et al.1 (|2006). We used results obtained by the MOPED 
algorithm, described in detail by Panteret al. (2007), to ex- 
tract star formation histories (SFH) from SDSS g alaxy spec- 
tra. It uses the stellar pop ulation models o f Bruzual & Chariot! 
( l2003l) and the|Chabrieil(|2003) initial mass function (IMF). 

The simulated galaxy catalog was obtained from the Mil- 
lennium Simulation (MS) of dark matter structure growth, 
which follows the evolution of = 2160'^ particles of 
mass 8.6 x 10^/i~^ M©, within a comoving box of size 
500fe~^ Mpc on a sid e, from redshift z — 127 to the present 
('Sprin gel et al.ll2005h . The adopted cosmological model is 
the concordance ACDM model with parameters = 0.25, 
fib = 0.045, h = 0.73, f^A = 0.75, n = 1 and erg = 0.9. 
These parameters are also adopted along this paper. Haloes 
and subhaloes in 64 output snapshot s were identified usin g the 
SUBFIND algorithm described in Springel et all ( 1200 Ik and 
merger trees were then constructed that describe how haloes 
grow as the model universe evolves. These merger trees form 
the basic input needed by the SAMs implemented in the sim- 
ulation. We used results obtained by the model described by 
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Fig. 1. — Mark correlation analysis for galaxies from the observational 
SDSS sample (filled circles; dotted lines) and from the Durham model (solid 
lines) and the Munich model (dashed lines). Error bars for the observational 
data were obtained with a jackknife procedure. There is an excellent agree- 
ment between the observed and the predicted correlations when both r-band 
luminosity (top panel) and stellar mass (bottom panel) are used as the marks. 



iDe Lucia & Blaizoi (12007 ) (Munich mo del), which is actu- 
ally a modified vers i on of m odels used in lCroton et alj (12006b 
and iDe Lucia et alJ (120061). and results f rom the model de- 
scribed in detail by iBower et alJ (l2006h (Durham model), 
which is an extension of t he GALFORM mod el implemented 
by lColeet aU (12000) and 'Benson et alJ (l2003l) . These SAMs 
have many differences in the overall procedure to form galax- 
ies from the dark matter halos. Particularly, the build up of 
merger trees, which form the basis of the SAMs integrated on 
numerical simulations, has been done in an independent way. 
Moreover, many aspects of the physics involved in the galaxy 
formation process are also treated in a different form, mainly 
those related to the inclusion of AGN feedback to regulate SF 
in galaxies. 

We derived the predicted galaxy mass assembly histories 
from the SAMs implemented on the MS. In order to compare 
with our nearby sample of SDSS galaxies, we selected all z = 
model galaxies with AB r-band absolute magnitudes Mr < 
—20, resulting in two base catalogs containing about 1.15 x 
10^ galaxies from the Durham model and about 1.81 x 10^ 
galaxies from the Munich model. For all galaxies in these 
catalogs, we recovered their complete merger trees by using 
the post-processing facilities in the MS databases. The mass 
assembly histories of the model galaxies were then retrieved 
by summing up the stellar masses of all the progenitors of a 
galaxy at each time-step of the simulation and the total mass 
assembled through SF between two time-steps was obtained 
by taking the difference of the total stellar mass retrieved at 
each time. Additionally, we have also computed the number 
of major mergers experienced by a galaxy along its entire life 
by summing up all merger events with mass ratio above 1/3 
of the main progenitor stellar mass at each time. 

3. MARK CORRELATIONS 

The SFH outputs, as well as other galaxy properties like 
luminosity and stellar mass, were used to compute the mark 
correlation functions for both the observational SDSS sam- 
ple and the MS galaxy catalogs. In our analysis, the mark 
correlations were computed with the help of the estimator 
WW/DD. Here DD refers to the number of galaxy pairs 
separated by distances r ± dr in each bin. While WW is the 
sum of the products of the galaxy properties weights W, also 
called marks, for the same pairs as in DD. The marks W 
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Fig. 2. — Mark correlations for 2 = galaxies from the catalogs generated 
by the SAMs but considering the positions of their main progenitors at 2 
2.4. It is clear from the bottom panel that galaxies that are massive today 
were more clustered in the past than less massive galaxies. 

are normalized so that the mean value (over the galaxies in 
the sample) is unity. This mark correlation ratio WW/ DD is 
therefore insensitive to the actual value of the mean cluster- 
ing (and its evolution). It is also insensitive to the mean value 
of the marks or its evolution. It just indicates if marks are 
spatially correlated above (or below) the average clustering. 

In Fig. [T]we show a comparison of the mark correlations in 
redshift-space for the SDSS sample and for the MS galaxy 
catalogs obtained by the models, where the marks are the 
r-band luminosity and the stellar mass. In the case of the 
SDSS results, the error bars in this and in subsequent fig- 
ures were estim ated by using a jackknife analysis, following 
IScranton et al.l (j2002), where we remeasured WW/ DD after 
removing a random region and repeating this procedure 30 
times. The comparison between the observed and predicted 
mark correlations is very good, confirming results obtained by 
previous studies which show that close pairs ten d to be more 
luminous and massive than the average (see e.g. ISkibba et al.l 
,2006.) . 

In the MS we can follow the evolution of galaxies by re- 
trieving the information stored in their merger trees. Assum- 
ing that the main progenitors represent the present-day galax- 
ies seen at different redshifts, we can investigate in which kind 
of environment a z = galaxy inhabited in the past. In Fig.|2] 
we do this exercise. We show the redshift-space mark correla- 
tions for z = galaxies in the model catalogs but considering 
their positions at z 2.4. Galaxy pairs containing massive 
galaxies today were more clustered, on average, at z ~ 2.4 
than their low-mass counterparts. The same trend is seen for 
luminosity, but to a less extent, for the Durham model and 
absent for the Munich model. Thus, in the simulations the 
relation between stellar mass and environment was already in 
place more than 10 Gyr ago. 

On the observational side, the recent advances in the field 
of spectral synthesis make it possible to retrieve detailed SFH 
from the spectra of nearby galaxies. Therefore, the evolu- 
tion of mass assembly in galaxies through SF as predicted 
by SAMs can be properly compared with results from obser- 
vational analysis. The mark correlation analysis can be then 
used to investigate the scale dependence of the SFH of galax- 
ies as predicted by the models and observations. In Fig. [3] 
we show results for both models and the SDSS data. In the 
top panels of this figure, the mark analysis were done con- 
sidering the total stellar mass formed (SFM) at a given look- 
back time bin as the marks. In the case of the SDSS, the 
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Fig. 3. — Top panels: Evolution of the weighted correlation functions as defined by the total stellar mass formed at different lookback time or redshift bins 
for the MS using the results from Durham model (solid lines) and Munich model (dashed lines), and for the observational results from the SDSS data using the 
MOPED (dotted lines; filled circles) approach to retrieve the SFH of galaxies. Bottom panels: MS merger analysis; the weights are associated to the number of 
major mergers (with mass ratio larger than 1/3) that the main progenitor of a z = galaxy experienced at each redshift. 



tr ends obtained in th is work reproduce the results discussed 
in ISheth et alJ (|2006). The most interesting result shown in 
Fig. [3] is the excellent agreement between the mark correla- 
tions obtained by the SAMs, especially the Durham model, 
and the observational data. The overall behaviour of the corre- 
lations in all lookback time bins shown in this figure indicates 
that close galaxy pairs as seen at z = formed more stellar 
mass at ^ 10 Gyr than the average, while more recently this 
trend is the opposite, with close pairs showing low levels of 
SF activity in comparison to the average. Note that for the 
first time bin, the Munich model predicts a positive correla- 
tion at smaller scales which is not caused by noise in the data, 
since the counts are similar to the Durham model. A possible 
explanation for these trends can be seen in the bottom panels 
of Fig. [3] where we show the mark correlations for the mod- 
els considering the number of major mergers experienced by 
a galaxy in a given redshift bin. There is a clear relation be- 
tween the increment in the SF activity of close z = galaxy 
pairs and the excess of major mergers at z > 1, thus reflect- 
ing that mergers played a significant role in the build up of 
galaxies in the high-z Universe. 

4. CONCLUSIONS 

In this Letter, we have investigated how the SF in galax- 
ies depends on the scale by comparing the mark correlation 
analysis of observed and model galaxies. Our results indi- 
cate that current semi- analytical gal a xy f ormation models, es- 
pecially the model by iBower et alJ (1200 6) (Durham model), 
implemented on the large Millennium Simulation are able to 
reproduce the trends seen in the mark correlations obtained 
for SDSS galaxies, when the marks are derived from the SFH 
of galaxies (as measured by the MOPED algorithm). Galaxies 
in close pairs seen nowadays formed stellar mass at a higher 
rate than the average population at z 2, and at lower red- 
shifts they show low levels of SF activity compared to the 
average. Using the predicted mark correlations from the mod- 
els, we find that for galaxies in z = close pairs the in- 
crease in their SFR at high-redshifts is intimately related to 
an excess in the number of major merger events they experi- 
enced in that epoch. From Fig. [3] it seems that the measured 
MOPED SF marks follow closely the marks of the numbers 
of mergers and seems to suggest that SF follows merger ac- 
tivity. It follows in this scenario that the observed downsiz- 
ing in SF properties of galaxies was driven by the environ- 
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ment primarily through mergers, which accelerated the evo- 
lution of massive galaxies seen today in close pairs (Fig. [TJ. 
The overall context of these results gives support to a natural 
path for galaxy evolution which proceeded via a nurture way, 
with denser environments playing a fundamenta l role in defin- 
ing ga lax y properties (see d etailed discussion in lMateus et al.l 
l2007l andlBundv et al.ll2006 , and references therein). 

The correlation properties of SF cannot be reproduced by 
models that only use the dark halo mass of the galaxy as the 
only parameter that drives the SF. This result is also related to 
the evidence of an environmental depe nden ce of halo forma- 
tion times found by lHarker et al.l (12006') and'Maul betsch et"an 
(2007). Actually, current implementations of the HOD 
paradigm do not include any environmental or spatial corre- 
lation. They can therefore not possibly account for the SF 
effect, which is clearly shown here to be an environmental ef- 
fect. It is therefore surprising that current SAMs (especially 
the Durham one) can reproduce the SF marks so remarkably 
well, given the complex processes that involve gas-baryon 
physics. Further, since we had no access to the simulations 
to fine tune their parameters, the obtained agreement is even 
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more remarkable. The strong correlation between the shape 
of the SF and the number of mergers mark correlations, which 
is our main finding, points toward a possible explanation for 
regulating the SF activity, which is something that galaxy for- 
mation models have inbuilt. However, if this was only the 
case the Munich models would be similar to the Durham ones. 
Their main difference resides in the way the merger process 
is translated to a SF recipe. In order to investigate this is- 
sue, in Fig. 2^ we show the evolution of the mean SFM in 
the same redshift bins of Fig. [3] where the scale dependence 
of the fluctuations relative to these mean values was previ- 
ously shown. Galaxies in the Munich model tend to form 
more stars at recent times than the Durham model but such 
difference is not too significative to produce the positive cor- 
relation shown in Fig. [3] Overall, the model galaxies form less 
stars at z > 0.5 compared to SDSS results and show an excess 
of SF in the first redshift bin. In Fig.|4j) we compare the dis- 
tributions of the recent SFM predicted by the models but now 
restricting to galaxies located in dense regions. These galaxies 
have local galaxy density values above the 25th percentile of 
the distribution of this quantity estimated through a /c-nearest 
neighbor approach (with k ~ 10). Galaxies from the Munich 
model which inhabit dense environments formed 3.3 times 
more stars in the first time bin shown in Fig. [3] than their 
counterparts from the Durham model. This higher recent SF 
activity predicted by the Munich model for galaxies in high- 
density regions could be the main factor responsible by the 
positive correlation at small scales shown in Fig. [3] This con- 
clusion is reinforced when we look the SFM distributions for 
the other time bins, which tend to be similar at high-z, where 



the mark correlations have almost the same values. Addition- 
ally, when we redo the mark analysis for the models consid- 
ering the same SFM distributions for galaxies in high-density 
regions, we obtain a lower correlation for the Munich model 
in the first time bin, confirming that the recent SF activity of 
galaxies in dense environments predicted by the models is a 
key factor which defines the differences in their mark correla- 
tions. This is also related to the way that models treat energy 
feedback from stars and AGN and how gas is made available 
to galaxies to cool and form stars (e.g Birnboim et al. 2007|; 
Dubois & Tey_ssjd 12008; Dekel & Birnboim 2008). Clearly, 
without further exploration of parameter space in the SAMs it 
is difficult to assert which parameter is the culprit of the good 
fit shown in Fig. [3] Current SAMs have no specific environ- 
mental dependence beyond a scale of 1 Mpc besides that of 
the merger process of dark matter and therefore their success 
must be linked to the way the merger process is then translated 
into gas forming stars. 



We gratefully thank the referee Darren Croton for helping to 
improve this paper This work was supported by the European 
Commission's ALFA-II programme through its funding of 
the Latin-american European Network for Astrophysics and 
Cosmology (LENAC). We warmly thank Andrew Benson for 
enlightening comments on semi-analytical galaxy formation 
models. The present work has been supported by grants from 
the Spanish Ministerio de Educacion y Ciencia, the Spanish 
CISC and the European Union FP7 program. 



REFERENCES 



Abazajian, K., et al. 2005, ApJ, 625, 613 

Angulo, R. E., Baugh, C. M., & Lacey, C. G. 2007. larXiv:0712.2280l 
Benson, A. J., Bower, R. G., Frenk, C. S., Lacey, C. G., Baugh, C. M., & 

Cole, S. 2003, ApJ, 599, 38 
Birnboim, Y., Dekel, A., & Neistein, E. 2007, MNRAS, 380, 339 
Bower, R. G., Benson, A. J., Malbon, R., Helly, J. C, Frenk, C. S., Baugh, 

C. M., Cole, S., & Lacey, C. G. 2006, MNRAS, 370, 645 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Bundy, K., et al. 2006, ApJ, 651, 120 
Chabrier, G. 2003, PASP, 115, 763 

Cid Fernandes, R., Mateus, A., Sodre, L., Stasifiska, G., Gomes, J. M. 2005, 
MNRAS, 358, 363 

Cid Fernandes, R., Asari, N. V., Sodre, L., Stasifiska, G., Mateus, A., Ton'es- 

Papaqui, J. P, Schoenell, W. 2007, MNRAS, 375, L16 
Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 319, 

168 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 

Croton, D. J., et al. 2006, MNRAS, 365, 1 1 

Dekel, A., & Birnboim, Y. 2008, MNRAS, 383, 119 

De Lucia, G., Springel, V., White, S. D. M., Croton, D., & Kauffmann, G. 
2006, MNRAS, 366, 499 



De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 
Dubois, Y, & Teyssier, R. 2008, A&A, 477, 79 
Gao, L., Springel, V., & White, S. D. M. 2005, MNRAS, 363, L66 
Marker, G., Cole, S., Helly, J., Frenk, C, & Jenkins, A. 2006, MNRAS, 367, 
1039 

Heavens, A., Panter, B., Jimenez, R., Dunlop, J. 2004, Nature, 428, 625 
Mateus, A., Sodre, L., Cid Fernandes, R., Stasihska, G. 2007, MNRAS, 374, 
1457 

Maulbetsch, C, Avila-Reese, V., Colm, P., Gottlober, S., Khalatyan, A., & 

Steinmetz, M. 2007, ApJ, 654, 53 
Panter, B., Jimenez, R., Heavens, A. R, Chariot, S. 2007, MNRAS, 378, 1550 
Scranton, R., et al. 2002, ApJ, 579, 48 

Sheth, R. K., Jimenez, R., Panter, B., Heavens, A. F 2006, ApJ, 650, L25 
Skibba, R., Sheth, R. K., Connolly, A. J., & Scranton, R. 2006, MNRAS, 369, 
68 

Springel, V., White, S. D. M., Tormen, G., Kauffmann, G. 2001, MNRAS, 
328, 726 

Springel, V., et al. 2005, Nature, 435, 629 



